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FOREWORD 


This  report  vas  prepared  at  the  Institute  of  Engineering  Research, 
University  of  California,  Berkeley,  under  Contract  AF  33(657) >8607, 

Project  7064,  "Aerotharaodynaaic  Investigations  In  High  Speed  Flow," 

Task  7064-01,  "Research  on  Hypersonic  Flow  Phenoasna."  Hm  work  was 
adalnlstered  under  the  direction  of  the  Aeronautical  Research  Laboratories, 
Office  of  Aerospace  Research,  USAF,  with  Capt.  Walter  W.  Wells  as  Task 
Scientist. 

The  theoretical  analysis  and  experlaental  work  were  carried  out  by 
H.  Tong  and  W.  H.  Giedt  under  the  supervision  of  Professors  G.  J.  Haslach 
and  S.  A.  Schaaf. 


11 


ABSTRACT 


Stagnation  point  haat  transfer  at  low  Reynolds  nuabars  in 
a  supersonic  air  streaa  was  investigated  ezperiaentally.  A  transient 
technique  was  eaployed  using  precooled  thin-walled  hemisphere-cylinder 
models.  Tests  were  conducted  at  nominal  Mach  nvnd>ers  of  2,  4,  and  6, 
and  in  the  Reynolds  number  range  of  80  to  1500.  Results  are  best 
represented  by  continuta  boundary  layer  theory.  Scatter  was  about 
ilO  per  cent,  so  that  any  possible  small  shock  wave-vorticity  effect 
at  the  lower  values  of  Reynolds  number  could  not  be  identified. 


TABLE  OF  COHTENTS 


I.  INTRODUCTION  1 

II.  THE  TRANSIENT  TECHNIQUE  4 

III.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE  6 

A.  Model  Construction  6 

B.  Instruaentstlon  8 

C.  Experiaentsl  Procedure  8 

IV.  EVALUATION  OP  HEAT  TRANSFER  COEFFICIENTS  10 

V.  DISCUSSION  13 

VI.  RESULTS  AND  CONCLUSIONS  16 

REFERENCES  18 

TABUS  I  20 

APPENDIX  22 


Iv 


LIST  OF  ILLDSTSATICIIS 


Figure  Page 

It  Coaperison  of  Theoretical  and  Bxperiaental  Stagnation 

Point  Heat  Transfer  38 

2.  Coaparison  of  Theoretical  Stagnation  Point  Heat  Transfer 

at  Low  Reynolds  Nuabers  39 

34.  Sectional  Drawing  of  Heat  Transfer  Model  40 

3B.  Instruaentation  Block  Dlagraa  40 

4  .  0.250  Inch  Diaawter  Heat  Transfer  Model  in  Test  Position  41 

5.  0.500  Inch  Diaaeter  Model  Before  Cooling  42 

6A.  Typical  TOaperature  Tiae  Record  43 

6B.  Typical  Seai-Log  Teaperature-Tiae  Plot  43 

7.  Coaparison  of  Measured  Stagnation  Point  Heat  Transfer  with 

Continuua  Boundary  Layer  Theory  44 

8.  Ratio  of  Average  Values  of  Measured  Stagnation  Point  Heat 

Transfer  to  Continuua  Boundary  layer  Theory  45 

9.  Typical  Hall  Teaperature  Distribution  at  Various  Tiaas  46 

10.  Coaparison  of  Mon-Disansional  Velocity  Gradients  as  Functions 

of  Mach  Muaber  47 


NOHENCIAIDSE 


c  specific  heat  at  constant  pressure 

P 

h  specific  enthalpy 

k  theraal  conductivity 

M  Hach  number 

Hu  Musselt  number 

p  pressure 

q  heat  flux  per  unit  area 

Pr  Prandtl  nuad>er 

R  gas  constant 

r  radius 

pud 

Ra  Reynolds  number 

T  absolute  temperature,  also 

U  velocity 

X  distance  from  stagnation  point 

constants 

5  s»del  wall  thickness 

7  specific  heat  ratio 

C  density  ratio  across  normal  shock 

X  constant 

|i  viscosity 

p  density 

T  time 


vl 


8ub«crlpt« 


o  atAgnation  reaarvolr 

1  upstream  of  shock 

2  doimstrssm  of  shock 

BL  boundary  layer  value 

e  edge  of  boimdary  layer 

FM  free  molecular  flow  value 

s  stagnation  behind  shock 

w  wall 


vll 


I.  INTRODUCTION 


The  flow  field  aroitnd  a  blunt  body  moving  through  the  atmosphere 

at  high  velocities  has  not  been  adequately  Interpreted  In  the  transition 

-“l  - 

region  between  the  continuum  and  free  molecular  flow  regimes  (1  >  ~  >  0.01). 

^1 

In  regard  to  heat  transfer  and  skin  friction  current  opinion  Is  that  a 
modified  boundary  layer  approach  may  be  possible.  If  second  order  effects, 
ordinarily  neglected  In  classical  boundary  layer,  are  considered.  In 
particular.  It  has  been  pointed  out  that  the  vortlclty  Introduced  as  the 
flow  passes  through  the  curved  bow  shock  wave  may  Increase  the  heat  trans¬ 
fer  and  skin  friction.  Other  second  order  effects  which  smy  enter  are 
slip  and  temperature  Jump,  boundary  layer  displacement  of  the  flow,  and 
body  curvature. 

Several  studies  of  this  transition  region  based  on  a  modified 

boundary  layer  analysis  have  been  reported.  Of  those  In  which  heat 

transfer  was  considered  two  distinguishable  approaches  were  followed. 

12  3 

The  first,  used  by  Hayes  and  Probsteln,  Hoshlzakl,  Ho  and  Probsteln, 

4 

and  Cheng,  was  based  on  treating  the  flow  between  the  body  and  the  bow 
shock  as  a  single  region.  Hayes  and  Probsteln  and  Hoshlzakl  assuate  the 
flow  between  the  shock  and  the  body  to  be  Incompressible.  Similar  results 
from  a  subsequent  study  by  Ho  and  Probsteln  Including  compressibility 
would  appear  to  justify  this.  Hoshlzakl  also  included  slip  and  temperature 
Jump  In  his  boundary  conditions,  but  found  their  effects  to  be  small. 

Results  are  presented  In  terms  of  stagnation  point  heat  transfer 
as  a  function  of  Reynolds  number.  To  assess  the  loiportance  of  the  second 
order  effects,  a  comparison  with  continuum  boundary  layer  theory  (References 
5  and  6)  Is  also  made.  Hayes  and  Probsteln,  Hoshlzakl,  and  Ho  and  Probsteln 
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obtain  similar  results,  which  Indicate  an  increase  over  continuum  boundary 
layer  (for  values  of  Re2  <  1000),  which  is  dependent  on  the  air  stream 
velocity  or  density  ratio  across  the  shock  (see  Figures  1  and  2).  Cheng, 
however,  predicts  an  increase  over  the  boundary  layer  theory  starting  at 
values  of  Re2  ■  20,000.  In  the  range  of  100  <  Re2  <  1000  his  results 
show  approximately  twice  the  increase  as  do  those  of  Hayes  and  Probstein. 
Also  of  interest  is  the  fact  that  after  reaching  a  maximum  his  results 
decrease  with  decreasing  Reynolds  nu^er  toward  free  isolecular  flow 
theory. 

The  second  approach  was  taken  by  Ferri,  Zakkay  and  Ting^  and  by 

0 

Van  Dyke.  Their  procedure  was  to  define  two  regions  between  the  shock 
and  the  body.  The  region  near  the  body  is  assuwd  to  be  viscous,  and  that 
near  the  shock  to  be  inviscid.  Boundary  conditions  include  specification 
of  switching  requirements  at  the  interface.  Both  analyses  predict  an 
increase  in  heat  transfer  over  continuum  boundary  layer  theory.  Van  Dyke's 
results  are  comparable  to  those  of  Ho  and  Probstein.  The  increase  predicted 
by  Ferri,  et.al.,  however,  is  from  200  to  300  percent  higher,  which  is  in 

general  agreeswnt  with  the  results  of  Cheng.  Stagnation  point  heat  transfer 

e 

data  which  agree  with  the  latter  theories  have  also  been  presented  by  Ferri, 
7  9 

et.al.,  *  .  This  was  obtained  by  using  a  transient  technique  in  a  blow-down 
type  tunnel.  ^ 

In  view  of  the  lack  of  agreement  between  the  theories  it  was 
felt  that  additional  experimental  results  were  needed.  The  low  density 
wind  tunnel  at  tha  University  of  California  can  provide  steady  flow 
conditions  in  the  upper  part  of  the  transition  region.  A  program  to 
obtain  beat  transfer  data  under  steady  state  conditions  was  therefore 
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conducted.  This  le  described  In  Reference  10.  The  results  obtained  are 


shown  In  Figure  1,  where  It  will  be  seen  that  they  agree  best  with  the 
theories  of  Van  Dyke  and  Hayes  and  Probsteln.  Also  Included  are  the 
recent  and  isore  applicable  results  of  Van  Dyke^^  for  a  free  stream 
Mach  nusril>er  of  4  and  a  wall-to-stagnatlon  temperature  ratio  of  0.5. 

Although  the  experimental  data  available  Indicated  substan- 
tlally  different  trends,  It  was  recognised  that  they  were  obtained  in 
different  test  facilities  and  using  different  techniques.  Because  of 
this  It  was  decided  that  some  additional  tests  In  the  low  pressure  wind 
tunnel  at  the  University  of  California  using  a  transient  technique  would 
be  desirable.  This  was  therefore  the  objective  of  the  present  study. 

The  basis  of  the  technique  Is  first  reviewed.  This  Is  followed  by  a 
description  of  the  isodels  developed,  the  testing  procedure  and  the 
results  obtained. 
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II.  THE  THAHSIENT  TECHNIQ^ 


In  gnncral,  this  method  for  determining  heat  rates  is  based  on 

e  knowledge  of  the  temperature  distributions  with  time  in  a  suitably 

designed  model.  Bvelustion  of  results  may  Involve  determination  of  both 

thermal  capacity  and  conduction  effects.  For  the  present  application  it 

is  convenient  to  utilise  a  thin~walled  model.  For  convective  heating  to 

a  thin-shell  hemisphere  the  temperature  distribution  is  governed  by  the 

12 

following  equation,  where  9  is  SMasured  from  the  stagnation  point. 

Here  the  first  term  represents  the  net  conduction  effect  and  the  second 
term  represents  the  rate  of  heat  storage.  Experimentally  it  is  desirable 
to  simplify  the  determination  of  the  heat  flux  by  evaltiating  q  irtien  the 
model  wall  is  isothermal.  Care  must  be  exercised  in  assuring  this  con* 
dltion,  as  it  will  be  shown  that  small  departures  from  the  isothermal 
condition  can  introduce  eubstantial  effects. 

Ta  use  the  transient  method  the  necessary  initial  temperature 
difference  between  model  and  free  streem  can  be  obtained  by  cooling  or 
heating  the  model,  raising  the  stagnation  tenperature  above  aiAient 
temperature,  or  a  coabination  of  these  methods.  In  this  experiment 
stagnation  temperature  heeting  was  not  possible,  so  that  the  choice  was 
between  heating  or  cooling  the  modal.  The  latter  alternative  was  used 
since  it  was  desirable  to  neglect  the  radiation  heat  exchange  which  for 
the  same  temperature  difference  would  be  greatest  for  heatingt 
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Itoo  procedures  were  considered  for  cooling  the  test  models.  The 
first  wes  to  fix  the  model  In  the  test  stresm  with  e  cooler  that  could  he 
rapidly  removed  from  the  air  stream,  or  alternatively,  the  model  could  he 
cooled  and  then  rapidly  moved  Into  the  air  stream.  Since  advantages  were 
not  obvious,  the  latter  alternative  was  used  because  of  the  simpler  mech¬ 
anism  which  could  be  employed. 

In  a  transient  heat  transfer  measurement  radiation  heat  exchange 
and  heat  losses  due  to  the  thermocouple  wires  must  be  considered.  The 
radiation  heat  exchange  can  be  made  negligible  by  using  a  polished  model 
and  small  temperature  differences.  The  effects  of  the  thermocouple  wire 
can  be  minimized  by  using  model  dlawnslons  such  that  the  effective  thermal 
capacity  of  the  model  Is  large  compared  to  that  of  the  wire.  This  can  be 
achieved  by  using  large  models  and/or  thick  walls  and  small  gauge  thermo¬ 
couple  wire. 
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III.  EXPERIMENIAL  EQUIPMENT  AND  PROCEDURE 
A.  MODEL  CONSTRUCTION 

The  nodalt  w«r«  foraad  by  •lectroplaclng  nickel  onto  an  alumlnua 
mandrel.  The  anode  material  wae  ordinary  cheat  nickel  bent  Into  the  ahape 
of  a  tube  to  that  the  mandrel  was  eurrounded  by  the  anode.  The  mandrel  was 
mounted  on  a  mechanism  which  imparted  a  reciprocating  smtlon  to  It  to 
achieve  more  uniform  deposits.  The  electroplating  bath  was  a  standard 
Watts  type  (240  gm  N1S0^»  45  gms  H1C1»  and  30  gms  Boric  Acid  per  liter 
of  water)  maintained  at  a  pH  factor  of  3  and  a  temperature  of  50*  C.  The 
plating  current  was  such  that  the  deposition  rate  was  approximately  0.001 
Inch  per  hour. 

The  models  were  easily  resovad  from  the  mandrels  by  cooling  the 
latter  In  cold  water.  They  were  then  mounted  on  a  Jig  and  two  0.004  Inch 
boles  were  drilled  at  the  stagnation  point  for  the  thermocouple  wires. 

40  gauge  copper •cons tantan  wires  were  passed  through  Individual  holes  and 
welded  In  place,  using  a  capacity  discharge  welder.  On  the  0.5  Inch  In 
diameter  smdel  for  which  local  temperature  measuresmnts  were  made  at  points 
ocher  than  the  stagnation  point.  Individual  copper  wires  were  welded  at 
various  locations  on  the  model  with  the  nickel  model  Itself  as  the  dls- 
slsdlar  metal.  A  second  copper-eons  tantan  pair  was  welded  at  the  base  of 
this  model.  Thus  It  was  possible  to  measure  the  tesiperature  of  the  baae 
and  the  temperature  difference  beMsen  the  base  and  any  desired  location 
on  the  model. 

Mo  satisfactory  data,  were  svallable  for  the  emf  of  copper-nickel 
Junctions;  hence  It  was  necessary  to  calibrate  the  model  Junctions.  This 
was  accomplished  by  mounting  the  smdel  on  a  liquid  nitrogen  cooled  collar. 
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which  was  in  turn  awuntad  Inslda  a  ball  Jar.  Iha  jar  was  than  avaeuatad 
to  a  praasure  of  about  40  Bicrons  of  narcurp  to  provida  a  controlled  tast 
condition.  Iha  stagnation  point  was  haatad  by  radiant  energy  from  a 
nlchroM  coll  to  provide  the  necessary  teaperature  difference  on  the 
tnodel;  thus  a  curve  was  obtained  for  the  copper-nickel  eaf  versus  the 
copper-constantan  eaf. 

The  aodels  were  aounted  on  fiberglass-reinforced  epoxy  stings 
which  had  approxlaately  the  aaas  coefficient  of  expansion  as  the  nickel 
aodels,  so  that  upon  the  cooling  the  aodel  to  atlng  joint  would  not  be 
under  excessive  stress. 

After  coapletlon  of  the  tests  the  aodels  were  cast  Into  an 
epoxy  resin  and  cut  In  half  along  the  axis.  This  surface  was  than 
polished  to  a  fine  finish  and  the  local  thickness  of  the  aodel  was 
aeasured  on  a  toolaaker's  alcroscope  with  a  stage  which  had  a  least 
count  of  0.0001  Inch.  Nlcro-photographs  were  also  taken  of  the  stag¬ 
nation  ragion  using  a  sMtallograph.  The  wall  thickness  could  thus  be 
deteralned  to  an  accuracy  of  1  percent. 

Three  aodel  sixes  were  used,  0.250  Inch,  0.500  Inch,  and  1.000 
Inch  dlaaeter,  with  0.004  Inch  noalnal  thicknesses.  Ttso  additional  0.250 
Inch  dlaaeter  aodels  were  forsied  with  wall  thicknesses  of  approxisately 
0.007  and  0.011  Inch  to  investigate  the  possible  effects  of  heat  losses 
due  to  the  presence  of  the  theraocouple  wire. 

A  sectional  drawing  of  a  sttdal  Is  shown  In  Figure  3A,  and  the 
0.500  Inch  and  0.250  Inch  disaster  aodels  as  sounted  In  the  wind  tunnel 
are  shown  In  Figures  4  and  5,  respectivaly. 
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B.  INSTRIMEMIATION 


Th*  thanocoupl*  mf  was  aapllflad  with  a  Waston  idllivolt-mlero* 
aoMtar  and  tha  variation  of  tanvaratura  with  raspact  to  tlaa  waa  racordad 
on  a  MlnnaapollS'Honaywall  oscillograph.  This  arrangsasnt  provided  a  com¬ 
bination  of  high  sensitivity  and  high  recording  spaed.  Sensitivities  as 
high  as  40  microvolts  per  Inch  were  possible  with  negligible  effects  from 
extraneous  noise  signals.  The  response  of  the  system  was  governed  by  the 
amplifier,  which  had  a  response  time  of  approxlmtely  0.1  second  for  a 
full  scale  step  Input. 

The  amplifier  used  Is  a  balancing  potentlomatrlc  Instrusmnt, 
and  thus  draws  no  current  from  the  signal  source.  This  simplified  the 
calibration  procedure  in  that  a  microvolt  potentiometer  could  be  used  as 
the  calibration  signal  with  no  correction  for  current  flow.  In  cases 
where  the  amplifier  was  not  used,  tha  calibrations  ware  performed  by 
ualng  an  auxiliary  voltage  source  In  series  with  an  equivalent  tharsm- 
couple  circuit  resistance  and  msasurlmg  the  smgnltude  of  the  source 
with  a  potentiometer. 

C.  EXPER1MBHT4L  PROCEDURB 

All  tests  were  conducted  In  the  Humber  4  Hind  Thnnel  at  tha 
University  of  California  Aeronautical  Sclancas  Laboratory.  This  tunnel 
Is  of  the  continuous  flow  type,  operated  by  two  five-stage  steam  ajactors 
with  a  msxlmisi  pumping  capacity  of  about  70  pounds  of  air  per  hour.  A 
complete  description  Is  given  In  Reference  13. 

Tha  test  models  were  mounted  on  a  high  speed  rotating  octopus 
such  that  tha  model  could  be  swung  Into  the  test  section  In  about  0.3 
sacond.  An  impact  pressure  probe  was  also  mounted  on  the  octopus  such 
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Mt  wta«n  the  model  tiee  out  of  the  test  section  the  probe  wss  In  the  test 
section.  In  this  way  pressure  measurements  could  be  made  while  the  isodel 
was  being  cooled. 

Cooling  of  the  models  was  accoogillshed  by  the  flow  of  liquid 
nitrogen  throt^h  a  copper  tube  wound  into  a  coil  of  approxlastely  1-1/4 
inch  Inside  diameter  (Figure  5).  One  end  of  the  coil  was  closed  and  the 
other  end  partially  closed  to  accomodate  the  a»del  to  be  cooled.  The 
cooling  coil  was  mounted  on  a  traverse  Mchanisn  which  permitted  it  to 
fully  enclose  the  model.  A  spring-loaded  trip  device  was  used  for  the 
rapid  removal  of  the  cooling  coil.  VRien  the  coil  was  fully  retracted, 
a  mercury  switch  was  activated,  which  started  the  rotation  of  the  octopua 
and  the  recorder  chart  drive.  To  speed  the  cooling  process  a  ssiall 
quantity  of  nitrogen  gas  was  bled  into  the  closed  end  ot  the  coll.  Thus 
between  tests  the  model  could  be  cooled  in  about  30  seconds.  A  typical 
experlsmntal  record  is  shown  in  Figure  6,  irtiere  a  chart  spaed  of  5  inches 
per  second  was  used. 
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IV.  EVALUATION  OF  HEAT  TRANSFER  COEFFICIENTS 

It  is  dsslxsbls  to  svalusts  tbs  boat  trsnsfsr  rstss  whan  tha 
■odal  la  Isotharaal  ao  that  tba  conduction  affects  ara  saro.  Whan  this 
is  tnia  the  taaparatura  variation  with  tlaa  will  be  axponantial;  there* 
fore  tha  data  ware  plotted  on  sewi'log  paper  and  tha  straight  line  portion 
was  used  to  evaluate  tha  heat  transfer  rata.  TO  confirw  that  this  corras* 
ponded  to  an  iaotharaal  wall  condition,  tanparatura  distributions  for  tha 
0.500  inch  diaaatar  wodal  ware  plotted,  k  typical  plot  is  shown  in 
Figure  6B.  It  was  found  that  tha  isotharaal  tlas  corresponded  to  about 
tha  aiddla  of  tha  axponantial  portion  of  tha  taaparatura-tlaa  plots. 

At  ths  tlaw  tha  nodal  was  isotharaal  tha  local  hsat  rata  could 
be  datamlnad  fron  tha  aquation 

-  -  P  «p  »  ^  (1) 

whara  tha  sadal  density  was  assunsd  constant  at  a  value  of  8.90  graas/ca  , 
and  its  specific  hast  as  a  function  of  tanparaturs  was  obtainsd  froa 
Rafaranca  14,  which  is  based  on  aany  published  data.  Tha  local  thickness 
was  obtained  fron  alcropbotographs  of  tha  nodal  cross-section,  as  described 
earlier. 

On  tha  saai-log  plots  ths  straight  line  portion  can  be  raprasantad 

w,  -  V  ■  •  V  • 

tuO 

where  O  is  tha  slope  of  tha  line.  Than 

S  -  -  V  -  ®  • 
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mnd  at  t  ■  0 


(I,  -  o 


t»0 


•ubsticutlng  Into  Equation  (1)  glvas 


(4) 


q  -  -  P  8  c  (T^  -  O 


T«0 


Foiming  tha  Nosaelt  nwbar 


Mu  - 


M 

k 


d  - a 

k  (i„  -  V 


t»0 


Hu  ■  -  ^  p  5  a 
K  p 


(5) 

(6) 

(7) 


Stagnation  point  thaoxy  baa  shown  that  tha  valoclty  gradlant  at  the  edge 
of  tha  boundary  layer  affects  tha  heat  transfer  rata,  and  It  Is  shown  In 
the  Appendix  that  this  valoclty  gradient  Is  dependant  only  on  the  Hach 
nuBd>er;  hence  to  account  for  the  gradient  the  Musselt  nui^er  Is  divided 
by  the  square  root  of  the  non-dlaenslonal  velocity  gradient 


1.681 


1  +  0.2 


(8) 


which  Is  obtained  (see  Appendix)  for  a  value  of  7  ■  7/5  and  a  pressure 
distribution  given  by 


P„  -  Pi  2 

.  1  -  cos*  0  .  (9) 

Pws  Pi 


11 


The  Reynolds  nusiber  Is  calculated  for  conditions  behind  the  shock 


and  is  given  by 


Re 


2 


P2U2 


d 


which  can  be  reduced  for  calculation  to 


(10) 


P2  d 

Re,  -  0.1685  ^  *  (11) 

1*2  ^ 


irtiere  P2  is  in  inches  of  vercury 

1*2  !•  in  1**“**«  (evaluated  at  ^) 

T  is  in  *R 
d  is  in  inches. 
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V.  DISCUSSION 


The  fMcasMxy  taaperature  dlffarcnea  b«tM«n  ths  test  ao4«l  wsll 
sod  tbs  stagnation  tsnparatura  of  the  tsst  air  straaa  was  obtained  by  cool* 
ing  the  iBodel  to  approxlnately  130*  K  prior  to  exposing  it  to  the  super> 
sonic  jet.  If  the  aodel  were  isotheraal  when  first  exposed  to  the  jet,  it 
would  be  sufficient  to  deterwine  the  gradient  of  teapereture  with  reepect 
to  tine  directly  froa  the  teapereture -tiae  records  by  saasuring  the  slope 
at  the  sudden  rise  in  the  record  which  represents  exposure  of  the  aodel 
to  the  jet.  Because  of  the  design  of  the  cooling  device,  however,  a  saall 
initial  teaperature  gradient  existed  between  the  90*  point  of  the  aodel 
and  the  stagnation  point.  The  errors  due  to  these  saall  gradients  have 
been  estiaated  (see  Appendix)  and  nay  be  significant,  it  was  therefore 
necessary  to  find  a  nethod  of  detemining  when  the  aodel  was  Isotheraal. 

Inspection  of  the  differential  equation  obtained  froa  a  heat 
balance  on  the  henispheric  portion  of  the  aodel  shows  that  when  the  aodel 
is  isotheraal  the  teapereture  is  an  exponential  function  of  tlae.  Re- 
plotting  the  teiBperature-tiae  behavior  on  a  seal-log  scale  revealed  that 
the  tenperature-tlae  behavior  was  exponential  for  a  short  period  of  tiaa, 
the  length  of  this  period  increasing  with  increasing  aodel  size  and 
decreasing  heat  rates.  Reinspection  of  the  differential  equation  indicated 
that  although  reaote,  it  was  possible  for  a  teaperature  gradient  to  be  such 
that  the  heat  transferred  by  conduction  could  also  be  exponential.  Thus  it 
was  not  necessarily  true  that  when  the  teaperature  was  exponential  with 
tine,  the  model  was  isotheraal. 

To  determine  the  teaperature  distribution  when  the  teaperature 
was  exponential  with  tine  the  0.500  inch  diaaeter  nodal  was  Instrumented 
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at  positions  othsr  than  tha  stagnation  point,  typical  distributions  at 
savaral  tisMS  ara  given  In  Figure  9,  from  vhlch  the  tlsM  idien  the  model 
was  Isotharmal  could  easily  be  determined*  this  tlaw  corresponded  very 
nearly  to  the  middle  of  the  exponential  part  of  the  stagnation  point 
temperature-tlsM  curves. 

T^ree  0.250  Inch  dlemeter  models  vlth  nominal  wall  thicknesses 
of  0.004  Inch,  0.007  Inch,  and  0.?11  Inch  were  used  to  chacfc  the  effects 
of  heat  losses  due  to  the  thermocouple  wires.  It  was  found  that  tha  heat 
transfer  rates  calculated  from  the  0.004  Inch  wall  model  were  approximately 
12  percent  less  than  those  calculated  from  the  0.007  Inch  wall  model  and  rates 
from  the  0.007  an&  0.011  Inch  wall  models  were  In  good  agreement.  A 
sesd-emplrlcal  analysis  was  used  to  estimate  the  errors  in  the  0.500  and 
1.000  Inch  diameter  models  (see  Appendix).  This  analysis  Indicated  that 
tha  errors  mould  be  small  enough  to  warrant  using  the  thin  walls  to 
achieve  larger  teaperature>tlme  gradients  which  could  be  s»re  accurately 
measured. 

On  tbs  basis  of  the  experience  gained  In  conducting  this  experl* 
ment  It  Is  felt  that  the  minimus  practical  slse  for  thin  wall  (0.004  inch) 
models  for  use  In  transient  heat  transfer  maasuransnts  Is  about  0.500  Inch 
diameter.  0.250  Inch  diameter  models  can  be  used  If  the  wall  thlcknese  Is 
approximately  doubled,  but  only  at  low  absolute  heat  transfer  rates.  It 
Is  also  essential  to  know  exactly  when  the  models  are  Isothermal  If  the 
calcttlstlon  Ignores  the  effect  of  conduction. 

For  future  Investigations  It  Is  recomsmnded  that  consideration  be 
given  to  locating  the  model  In  the  stream  and  suddenly  removing  a  cooling 
shield.  It  may  be  possible  In  this  way  to  achieve  a  faster  exposure  of  the 
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nodel  to  the  alritreea.  Attention  should  be  given  in  the  design  of  the 
cooling  shield  so  that  it  will  aaka  the  aodel  isotheraal.  Also  in  this 
regard  it  is  noted  that  heating  of  air  flowing  through  the  tunnel  would 
provide  additional  advantages.  Higher  initial  tesiperature  differences 
would  be  possible,  and  aodels  could  be  cooled  by  allowing  thea  to  cosw 
to  an  isotheraal  condition  at  rooa  tes^erature.  Developaant  of  this 
latter  suggestion  should  iaprovs  the  accuracy  with  which  the  initial 
slope  of  the  aodel  tesiperature  change  can  be  deterained. 
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VI.  RESULTS  AND  CONCLUSIONS 


The  stagnecion  point  heat  transfer  reaults  determined  as  described 
in  previous  sections  are  listed  in  Table  I.  They  are  shown  graphically  in 
Figure  7  in  terms  of  the  Nusselt  number  divided  by  the  square  root  of  the 
dimensionless  velocity  gradient  at  the  edge  of  the  boundary  layer,  Nu/  >/3  *, 
as  a  function  of  the  Reynolds  number  downstream  of  the  shock,  Re2> 

Generally  several  runs  were  made  at  each  test  condition.  Evaluation  of 
all  the  records  thus  obtained  yielded  duplicate  results  for  most  values 
of  Re2.  Since  the  variation  of  Nu/  obtained  at  any  value  of  Re2 

was  indicative  of  the  experimental  accuracy,  all  data  points  were  Included 
in  Figure  7.  As  can  be  noted,  only  a  few  points  fell  outside  the  lines 
deviating  llO  per  cent  from  the  continuum  boundary  layer  reference  line. 

An  average  Nu/  >/u^'  was  computed  at  each  value  of  Re2  for 
which  two  or  more  experimental  points  were  available.  The  ratios  of  these 


average  values  to  continuum  boundary  layer  theory  have  been  plotted  as  a 
function  of  Re2  in  Figure  8.  The  curve  representing  a  least  squares  fit 
of  the  data  reported  in  Reference  10  has  also  been  included  in  this  figure. 
Although  the  results  from  this  Investigation  are  In  agreement  with  this 
curve.  It  should  be  noted  that  the  data  at  the  lower  Reynolds  nua8>ers  are 
priamrlly  from  the  0.250  Inch  sndel  tested  at  ~  4.0.  In  view  of  this 
and  the  experiiaental  scatter,  a  firm  conclusion  regarding  an  Increase 


over  boundary  layer  theory  is  not  considered  Justified. 

The  results  are  therefore  considered  to  be  In  agreesient  with 
boundary  layer  theory  over  the  range  of  Reynolds  numbers  Investigated. 
However,  since  the  trends  predicted  by  the  analyses  of  References  1,  2, 
3,  and  8  lie  within  the  experlatental  scatter.  It  must  be  recognised  that 
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small  second  order  effects  may  be  present.  Reference  to  Figure  1  shows 

7  9 

that  the  experimental  results  of  Ferrl,  et.al.,  '  as  well  as  the  theory 
7  4 

of  Ferrl  and  Cheng  are  significantly  higher,  and  also  show  a  deviation 
from  continuum  boundary  layer  theory  at  a  much  higher  Reynolds  number. 

It  Is  noted  that  the  data  of  Ferrl,  et.al.,  were  obtained  at  higher  stag* 
nation  temperatures  (**  2300  R  compared  to  530  R).  This  has  been  advanced 
as  the  reason  for  the  differences  in  the  theoretical  predictions.^^  This 
explanation  Is,  however,  questioned  by  Van  Dyke.^^ 

In  reviewing  the  findings  of  this  Investigation  it  is  concluded 
that  the  transient  technique  developed  for  stagnation  point  heat  transfer 
measurements  yields  results  comparable  In  accuracy  to  the  steady  state 
technique  used  previously.  Data  obtained  are  In  general  agreement  with 
the  previous  results^^  throughout  the  available  Hach  and  Reynolds  number 
ranges.  In  general,  the  data  scatters  about  110  per  cent  above  and  below 
continuum  boundary  layer  theory,  but  tend  to  cluster  above  at  the  lower 
values  of  Reynolds  number.  Although  the  predictions  of  References  1,  2, 
3,  and  8  are  Included  In  this  scatter,  no  definite  conclusion  regarding 
second  order  effects  can  be  reached. 
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TABLE  I,  TABULATED  RESULTS 


1  inch  Hodel 


Re2 

Nu 

Nu/ 

^/‘^BL 

Run  A 

Run  B 

Run  C 

Run  D 

6.01 

1485 

62.0 

25.6 

27.1 

26.7 

0.993 

5.98 

1260 

56.2 

23.2 

24.7 

0.980 

6 

5.90 

1071 

51.0 

18.6 

23.2 

22.0 

0.944 

o 

5.82 

866 

44.9 

18.6 

19.8 

0.944 

5.64 

640 

44.0 

18.3 

16.6 

1.080 

4.00 

1 

41.5 

■ 

■ 

1.000 

3.98 

39.1 

1.017 

~  4 

3.92 

520 

37.5 

^^9 

15.8 

mM 

3.86 

418 

31.1 

13.5 

15.0 

|M 

0.985 

3.70 

314 

27.0 

11.9 

IB 

1.001 

2.22 

1019 

44.5 

22.3 

2.20 

681 

36.8 

18.5 

1.058 

-  2 

2.18 

549 

31.8 

16.1 

1.028 

B 

2.18 

549 

32.9 

16.7 

1.065 

2.15 

369 

27.1 

13.8 

1.070 
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TABLE  1  -  Continued 


1/2  Inch  Model 


APPENDIX 


BOUNDARY  LAYER  HEAT  IRANSFER  REUTION 

Proa  Fay  and  Rlddall^  at  tha  atagnatlon  point 


Mu 


<1 

1  Pr 

1 

1 - 1 

m 

m 

- 

1/2 

0^ 

• 

tl 

w 

• 

'^nd 


"«x  -  0.67  I 


P  P 

wt  ws 


0.4 


where  Mu 


q  X  r 


X  k(h  'h  ) 
'  o  wa' 


and  X  la  tha  local  coordinate. 

Equating  Equationa  (A>1)  and  (A-2)  givaa 

0.4 


0.67 


P  P 
o  o 


172 

1 

^wa  **Wa  (  dx  }  , 

m  * 

for  low  atagnatlon  taaperaturaa  pp*-  conatant  ao  that 
and  replacing  P^,P^,  by  P2P2 


0.67  - 


<1  _ 

r 

ri*!  (s^). 


for  a  perfect  gaa  at  the  atagnatlon  point 
<*‘o-*W>  -  «p  <^o"^> 
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(A-1) 


(A-2) 


••  1 


(A-3) 


and  since  Pr  ■ 


®t»2 

^  ^  t  Equation  (A-3)  can  be  written 


0.67  - 


(VV^ 


^2 


lV2 


2  /‘*“e 
dx 


Defining  the  Nuaselt  and  Reynolds  nund>ers  as 


-  h  d  q  d 

r”  "  7t  .t  r  k 

*2  ^ws^  “2 


Ra,  ■ 


P2U2d 


we  obtain 


0.67  - 


Nu 


or 


Mu 


.'d  (x/d) 


17? 


du  \ 

d 

0- 

*1  I 

2 

'  s 

17? 


•  0.67  ■'/ Re- 


d  (u^/Uj) 

Defining  "j  (x/d)  "  '*e'  *  above  expression  bec< 


-  0.67  (A-4) 

e 

in  which  all  properties  and  the  velocity  are  based  on  conditions  behind 
the  normal  part  of  the  shock.  Equation  (A>A)  is  in  agreeoient  with  the 
incoBipressible  theory  of  Slbulkin.^ 
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D1MENSI(H(LESS  VELOCITY  GRADIENT 


From  Newton's  Second  Law 


du 

u  — £  . 

1 

‘^Pe 

dx 

dx 

Differentiating  once  with  respect 

to  the  X 

coordinate 

d^u 

/du 

e 

1  2 

L.  ^ 

dp^ 

“e  ^  2  +  1 

dx 

[dx  ^ 

P  .2 

e  dx 

+ 

dx 

(B-1) 


At  the  stagnation  point  u  ■  0  and  p  «  constant  ■  P  ,  and  due  to 

6  6  6 

syiMetry  (dp^/dx)  ■  0  ,  so  that 


du 
_ c 

dx 


(B-2) 


For  a  henlsphere  the  pressure  distribution  is  given  by  the  modified 
Newtonian  relation 


hUi 

Ps  ■  Pi 


2 

cos^  e 


since  by  the  boundary  layer  concept  p^  **  p^ 
Differentiating  and  noting  that  x  ■  rG 


2 

“  “  (P,  -  pj)  cos  0  sin  0 

-  ^  -  (p^  *  pj)  [cos^  0  -  sln^  0] 

at  the  stagnation  point 

cos  0  ■  1  end  sin  0  >  0 


dx 

-V 


dx‘ 
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so  that 


s 


substituting  into  (B-2)  yields 


w  I  9 

jet)  ,  ■  77 -  "!> 


From  tha  previous  section 


.  du 

U  •  -  ~  T-^ 


V  ■  S;  I 


2  r  2  / 

-  Bj  [  P,  <>“. 


for  a  perfect  gas  ^  »  RT 


Z  •  .  2  ■/a 

“e  “  Oj 


RT  -  RT  -r 

L  •  •  Ps 


r  Pj 

«2  47  L  "Ps 


2*f2  .if,  li  1  h$ 
"  ”2  V7L  PsJ  Jh 


(B-3) 


assuming  ~  T2«  Substituting  the  relations 


J±L 


27Mj  -  (7-1) 


and 


2  1 
Mj  (7-1)  +  I  ; 

L  27Mj^  -  (7-1 )j 


into  E<iuaClon  (B-3)  leads  to 

S'  - 

e  V  7 


2711^^  -  (7-  1)]'/' 


Mi(7-1)  +  2 


7-H 


27Mj  -  (7-1) 


1/2 


2  -|l/2  j-  27M,^  -  (7-1)  -  (7+1)''^^^ 


27Mj  -  (7 

Lm^Vi)  +  2  j 


27Mj‘  -  (7-1) 


-  2^2 


M, 


1  I 


1/2 


1  +  ^ 


rVJ 


Finally  for  7  ■  1.4 


2.88 


Mi^-l 


1  -f  0.2  M 


1  J 


1/2 


to-4) 


Probstein^^  defines  the  non-diaaensional  velocity  gradient  as 


which  can  be  written  as 


••  • 
u 

a 


CB-5) 
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Probtteln  and  Kenp^^  define  the  non-dlawnslonal  velocity  gradient  as 


which  can  be  written  as 


Relations  (B-4),  (B>5),  and  (B-6)  are  plotted  In  Figure  10  for  coaparlson. 
Equation  (B*4)  was  used  In  the  calculations  for  this  ezperlaent* 
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ESTIMAIB  OF  HEAT  FLUX  ERROR  DUE  TO  SLOPE  MEASURBfEHTS 


A  simple  fin  solution  for  s  hsmlsphers  csp  yields 


q  -  + 


h  6  r 

I 


cos  9 
sin  6 


-  »  ®p  ^ 


irtisre  T  ••  stagnation  point  this  reduces  to 


q 


+ 


2  k  6  d^T 
r2  ae" 


p  6  Cp 


ar 


<c-l) 


For  the  conditions 


k  -  0.7  5*^^=  ,  5  -  0.010  cm  ,  pB  -  0.089  ®V  * 

cm  K  cm 

2  2 

c  *0.3  36ules/gm  *K  »  and  r  •  0.403  cm  <0.5  inch  model). 


’  ■  +  "  ’‘o?w  ^ 


a^T 


-  +  0.348  ^  -  0.0267 

ae* 


iC-2) 


For  the  0.500  inch  model  at  Re  *  320  and  M  -  6  the  exponential 
portion  of  the  tesiperature>tims  curve  occurs  at  t  ••  0.05  sec.  This  has 
been  verified  to  be  the  time  when  the  smdel  is  isothermal.  At  T  ■  0.10 
sec  the  temperature  difference  between  the  stagnation  point  and  the  40* 
point  is  approxlmatsly  0.3  **.  If  we  assume  a  cosine  temperature  distri¬ 
bution 

^  '  V.'  ■  c.  S  <C-J) 

^  •  ’w* 


28 


-  ’so*>  - 


0.3  -  (I  -  0.767)<^  -  T^,) 

\n-h0‘  ■  *•*“ 

Thar* for* 

®  1 .  ■ 

and  for  th*  statad  condition*  *  *  27  *l^*ae. 

Subatltutlng  Into  (C-2)  glvaa 

q  -  1.288  X  0.0348  -f  0.0267  x  27 

■  (0.045  +  0.721)  watt/oi^ 

Fxob  this  calculation  «a  obsarva  that  th*  conduction  tara  can  ba  a 
larga  (4.8  parcant  In  thla  caaa)  portion  of  th*  total  haat  flua  for 
ralatlvaly  aaall  taaparatura  gradlanta  along  tha  aodol  aurfaca. 

In  th*  raductlon  of  tha  data  tha  taaparatura- tlM  alopaa  oar* 
avaluatad  at  th*  atralght  llna  portion  of  th*  aaal-log  plota  and  th* 
larga  arrora  aaaoclatad  vlth  aaall  dovlatlona  froa  thla  tlaa  account 
for  aoaa  of  th*  acattar  In  tha  azpariaantal  raaulta.  It  la  aatlaatad 
that  th*  arrora  du*  to  alopa  aaaaurMMnta  ara  of  th*  ordar  of  5  -  6 
parcant. 
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BSTDttkn  or  BADUTIOI  HEAT  BXCHAHBB 


▲••OBliit  that  tte  t«sc  aodal  la  gray  and  tha  ehaabar  a  black 
iaothazaal  aneloanra,  tha  radiation  boat  osehanga  can  ba  axpraaaad  aa 

^  ^  4  _  ,  ,,  ,a“W  aatta 

nhara  the  aubacrlpta  1  and  2  rafar  to  tha  body  and  tha  aneloaura, 
raapaetlvaly.  For  pollabad  nlekal  tba  aalaalvlty  for  nomal  radiation 
la  approxlaataly  0.045  -  0.087  (aao,  for  Inatanca,  Harka*  Handbook) 
and  for  hanlapharle  radiation  thla  valua  la  Ineraaaad  about  20  poreant 
ao  that  an  anlaalvlty  of  0.1  would  fom  an  tvpar  bound  on  tha  haat 
aaebangad  by  radiation.  For  anblant  and  wall  taaparaturaa  of  300*  K 
and  130*  JL,  raapaetlvaly 

q  -  0.10  X  5.71  X  10‘“  [300^  -  130^) 

*  0.0044  watta/en^ 

2 

Tha  lowaat  beat  rataa  naaaurad  vara  approKiaataly  0.4  watta/cn  •  ao 
that  tha  naataiM  corractlon  la  of  tha  ordar  of  1  pareant. 
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BRBOR  DUB  TO  PSB8ERCB  OP  THBSHOCOUPLB  WIRB 


A  hMt  b«Une«  on  an  alaaant  of  a  ehanaocoupla  wlra  ylalda 


‘QZZD — ^  "*2 


’1  •  -  ».  I*  If  \  ’2  •  -  i  >j 


’l  -  ’2  *  ’3  -  "  'p  \  ^  ^ 


Thaxafora 


M  >,  ^  P  A* 


or 


U- 


AsauadLng  a  aoluclon  of  tha  form  and  lattlng 


T  -  X(x)  0  (T) 


9m 


glvaa 


S!1 

X 


2  0*  2 
-  P  ^  ■  constant  •  \ 


IharafoTOp 


dn  0  "  “  “T  0  +  dn  C. 

r  ^ 


or 


0  ■  Cj  e 
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and 


X"  -  X‘  X  -  0 


The  general  aolution  of  this  equation  la 


Thus 


+  Xx  -  Xx 

X  ■  C2  a  +  Cj  e 


^aw-*^ 


«  e 


Xx  -  Xx 

(C2  a  +  Cj  a  ) 


dr 

and  ^  «  a  XCCj 


Xx  -  Xx 

a  -  C3  c  ) 


Since 


Also 


^—0  Cj.O 


dr 

3F 


p2  «3  • 


At  X  •  0 


x^ 


-  Cj  X  e 


and 


-  -  C,  —  a 

7tf^  ^  fT 


For  the  highest  heat  transfer  conditions  ^  *  0.25  and  (  ^  ) 

P*  ^  t-0 

•  34.4  *K/sse,  bancs 
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=3  ■  srn 


»  135 

(  ^  )  ^  -  135  •%/cm  . 

T*0 


Now  for  0.003  inch  wires 


-  ‘JxpO  -  -  2  k  f 
T»0 


X  <  ^ 


2  X  0.7  2a£|-  X  4.55  X  10“®  cn^  x  135  — 

CWk  K  CB 


>  0.0086  watt  . 


If  we  consider  e  disk  with  an  area  equivalent  to  an  arc  of  the  hasdsphare, 
the  heat  flux  stored  In  this  aleaent  is 


Q  <■  p  c  &  A  ^ 
2  p  s  or 


«.  p  c  B  X  (r  sin  0)^  ^ 


Q  ■  Qi  +  Qa 


Vi 

The  error  is  approximstely  ^  and  for  the  1/4  inch  aodel  this  was  found 


to  be  *>  12  percent.  Estiaating  Q2  froa  the  aeasured  value  of  dl/dr  gives 


Q2  -  8.9x0.3x0.01x3.14x0.101  sin*  0  x  34.4 


-  0.292  sin  0 
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Hm  value  of  0  It  now  dttemlntd  by  noting 


Qi 

Q1  +  Q2 


0.12 


and 

Qj  -  0.12  (Qj  +  Qj) 

<*2  ■  0?S  -«l  - 

Therefore 

0.292  tln^  0  -  7.34  x  0.0086 
tin^  0  -  0.216 


tin  0  ■  0.465 

0  -  26.8* 

Ataualng  that  the  taae  turface  area  la  Involved  (and  wall  to  wire 
thlckneaa  which  wat  held  conatant),  for  the  1/2  Inch  andel  for  the 
taae  condltiona 


Qj  "  0.544  watt 

and 

2l  0.0086  . 

Q  "  0.0086  +  0.544  * 

Siailarly  for  the  1  Inch  aodel  the  error  la  found  to  be  ■'  0.4  percent. 

From  thla  approxlatte  calculation  we  conclude  that  the  reaulta 
at  preaentad  aay  be  ayateaatlcally  low  by  1.7  and  0.4  percent  for  the 
1/2  Inch  and  1  Inch  dlaaeter  aodela,  raapectlvely. 
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ERROR  ANALYSIS  FOR  INSTRUMENIATION 


The  heat  trenafer  coefficient  Is 


Nu  - 


k  P  6  Cp  « 


Ihe  errors  due  to  the  quantity  dp  6  are  systematic  and  the  values 

18 

of  c^  and  k  are  obtained  from  published  data  and  assumed  accurate. 
The  remaining  quantity  is  a  so  that  the  maximum  probable  error  in  this 
value  is  the  maximum  probable  error  In  the  Nusselt  number.  From  Section 
IV  we  have 


dT 


T»0 


-  ■'t.o’  “ 


dt 

dx 


T-0 


«o  -  Vo* 


The  probable  error  can  be  written  as 


ot 

a 


®(dT/dT) 

2 

dt/dT 

X 

W'W  J 


n  1/2 


The  errors  in  the  slope  measurements  have  already  been  estimated  to  be 
~  6  percent.  The  error  in  the  temperature  measurements  depends  on  the 
calibration  source  and  tbs  ability  to  read  the  records.  This  error  is 
probably  less  than  3  percent.  Hence 


^  -  [0.0036  +  0.0009]^/^ 

-  6.7%  . 
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From  Section  IV  the  Reynolds  number  Is 


P2  Mj  d 

Re^  -  0.1685  —  — 

The  errors  in  M2  are  negligible  since  behind  the  shock  the  Mach  number 
is  a  very  small  function  of  the  pressure  ratio  p.  /p  .  The  error  in 

O 

d  will  be  systematic,  so  that 


Re2  -  K  — 
^^2 


where  K  is  a  constant.  Assuming  a  linear  viscosity  law  **  T2  and 
noting  that  at  a  given  Mach  number  T2  *>  we  have 


Re. 


K 


«o> 


3/1 


The  probable  error  can  be  written  as 


All  iaq>act  pressures  were  measured  with  an  oil  aanosMter  with  a  least 
count  of  0.001  inch.  Stream  variations  with  tlsm  are  about  twice  this 
smgnitude.  The  lowest  impact  pressures  smasured  were  **  0.500  inch,  so 
that  the  error  due  to  this  source  would  be  0.4  percent.  Stagnation 
temperature  smasurements  could  be  read  to  within  1.0*F,  and  for  an 
ambient  temperature  of  540*R  the  error  is  *  0.2  percent.  Hence 
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bs. 

Re 


[8 


10“®  4.50  +  10 


-  0.5X  . 
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F16DRB  1.  OOMPARISOII  OF  THEORETICAL  AMD  EXPERIMENTAL  STAOIATIOM  POINT  HEAT  TRANSFER 


FI6URB  2.  CGHPARISON  OF  IHBORXnCAL  STMaUTZON  POINT  HB4T  TRANSFER  AT  LON  RBTNOIOS 


Figure  3B.  InatruMiiteclen  Block  DUgraa 
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FIGURE  4.  0.250  INCH  DIAMETER  HEAT  TRANSFER  MODEL  IN  TEST  POSITION 
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FIGURE  5.  0.500  INCH  DIAMETER  MODEL  BEFORE  COOLING 
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FICDSE  9.  TYPICAL  WALL  TEMPERATURE  DISTRIBUTIOM  AT  VARIOUS  TIMES 
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OP  mH-OnBISIOlUL  VILOCm  GBADXBNTS  as  FOMCnOMS  OF  MACH  MOMBBR 


I 

i 

I 

I 
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